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The letter resolves several intriguing and fundamental aspects of debonding at soft interfaces, 
including the formation and persistence of regularly arranged nanocavities and bridges, "adhesion- 
debonding hysteresis", and vastly lower adhesive strengths compared to the absence of pattern 
formation. The analysis shows the hysteresis to be caused by an energy barrier that separates the 
metastable patterned configurations during withdrawal, and the debonded state. The metastable 
morphological pathways involving cavitation and peeling of contact zones engender substantially 
lower debonding forces. 
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The related phenomena of adhesion, debonding and 
interfacial cavitation or cracking at soft elastic inter- 
faces have been intensely studied both in view of their 
technological applications and the many unresolved sci- 
entific issues related to the pathways, morphology and 
forces of debonding. A rigid surface (contactor) ini- 
tially in contact with a soft elastic film, upon with- 
drawal, debonds by the formation of a pattern of well 
defined spacing consisting of areas of intimate contact 
and interfacial cavities 0^0, 0, H- A linear stability 
analysis!^ showed spontaneous surface roughening of 
the film when the contactor is initially brought in close 
proximity (<20nm). The wavelength (A) of this sur- 
face pattern depends only on the film thickness {h) 
(A ~ 3/i), but is independent of the strength and na- 
ture of the adhesive interactions as well as the elastic 
properties of the film. The regions of adhesive contact 
persist until the contactor is pulled to a much larger 
distance than that at which the initial contact insta- 
bility was triggered upon approach. This difference 
between approach and withdrawal behaviors may be 
referred to as "adhesion-debonding or contact hystere- 
sis". Another interesting aspect is the pull-off force 
required for debonding can be several orders smaller 
than the force calculated based on the assumption of 
flat surfaces. Clearly, such a significant reduction can- 
not be explained merely by ~50% reduction in the con- 
tact area observed at detachment 0- The mechanism 
of "adhesion-debonding hysteresis" , its associated mor- 
phologies including patterned cavitation, pull-off force 
and distance are the key unresolved issues in debonding 
at soft interfaces that are addressed here. 

Figure. ^ illustrates the film-contactor geometry. 
The total energy consists of the stabilizing stored elas- 
tic energy and the destabilizing attractive interaction 
between the contactor and the film. 



(1) 



Contactor (Approach) 



Contactor t (Withdrawal) 



n= / W{e)dV+ / U{dQ-u-n)dS. 



V 



s 



W is the strain energy density defined as, 
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FIG. 1: Left: Schematic of small amplitude pattern forma- 
tion on approach at critical distance dc- Right: Column 
stretching and growth of cavities on withdrawal of the con- 
tactor. 



(<10MPa). The interaction potential consists of an at- 
tractive van der Waals component along with a short 
range Born repulsion, represented by U (do — u ■ n) = 
-A/127r(do - u ■ + B/(rfo - u ■ n)^ where A is 
the Hamakcr constant (of the order of 10~^°J) and 
B is the coefficient of Born repulsion. The coeffi- 
cient B is correlated to the adhesive energy at contact 
{AG = U{de) = A/{16TTdl)), where 4 is the equilib- 
rium separation distance obtained from U'{de) = 0. 
This form of interaction implies that the force required 
to pull off two rigid flat surfaces is F^l^^j. = AG/de- 
As may be expected, our detailed studies (to be pub- 
lished) have confirmed that the debonding is controlled 
by the adhesive energy rather than the detailed func- 
tional form of the potential. Based on the linear sta- 
bility of eqn. it was shown that the film sur- 
face becomes spontaneously rough as the contactor ap- 
proaches it to within a small critical distance dc at 
which h\U"\/fi > 6.22. The lengthscale of the pattern 
consisting of cavities and bridges (regions of contact) 
is about three times the film thickness regardless of the 
interaction potential and the elastic properties both of 
which are in agreement with observations [2 lilJI- 

We explain the physical basis of the adhesion- 
debonding hysteresis (debonding distance 3> dc) by 
the following simple analysis which is also in confor- 
mity with the detailed simulations presented here. For 
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FIG. 2: Energy landscape as a function of wavenumber. 
The plots show the presence of an energy barrier separating 
the flat film and patterned state for large values of do/de- 
The physical parameters are /i = O.lMPa, h — lOmicrons, 
A = 10-2° J, AG = l.OmJ/m^ 



a single Fourier mode, U2{x) — coskx, the total en- 
ergy (per unit length of the film) is shown in fig. |21 
for two different values of the gap-thickness do above 
and below the critical distance dc- As predicted by 
the linear stability, the patterned configuration with 
hk « 2.12 (A ~ 3h) has the lowest energy for do < dc 
rather than fiat film. However, for do > dc, the fiat 
film configuration has the lower energy, although the 
patterned state remains a local minimum, metastable 
state. For do > dc, the global (fiat film) and the lo- 
cal (patterned state) minima are separated by a large 
energy barrier (fig. It is due to the presence this en- 
ergy barrier that the patterned states formed during the 
approach persists in its metastable configuration upon 
withdrawal. Of course, as shown by simulations, the 
pattern is far more complex in that it consists of many 
Fourier modes, leading to a multiplicity of metastable 
states of varying energies. Thus, during the pull off, 
the system "hops" through a succession of metastable 
states leading to a strong "path dependence" especially 
in the presence of heterogeneities, noise etc. 

The complete simulation of the pull-off process is 
achieved by the Fourier representation U2{xi,0) = 
^^Jg"^ a„ cos(fc„a;i), where a„ is the amplitude of the 
n-th Fourier mode with wavenumber fc„(= 2'Kn/L). 
The total energy per unit depth of the film 
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The stresses that develop in the film are 

determined from the Fourier coefficients as 

0-22(3^1,0) = 2nJ2n=0 "-nkn S{knh) cos knX I, whcrC 

S{^) = ^-^-^^^^z^-- F is the average force per unit 
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FIG. 3: (A) Film profiles at various separation distances 
without external noise and step size of Inm [A = IQ-^^J, 
/i = O.SMPa, h = lOmicrons) (B) & (C) Variation, re- 
spectively, of fractional contact area a and force F/ fi with 
separation distances: Curves (1), (2) and (3) correspond to 
the step size of O.lnm, Inm and 2nm, respectively. Curves 
(1) and (3) correspond to the catastrophic column collapse 
and continuous peeling modes of failure, and curve (2) is an 
intermediate. (D) Normal stress distribution along the sur- 
face of the film for cases in (A) showing maximum stresses, 
responsible for peeling, at the column edges. 



area exerted on the contactor plate to hold the film in 
equilibrium at a given separation. 

A conjugate gradient(CG) scheme (which finds the 
local minimum closest to the initial configuration) was 
employed to find the Fourier coefficients that result in 
a minimum energy pattern for a given separation dis- 
tance do > dc starting from the contact proximity. The 
robustness of the energy minimum thus isolated was 
confirmed by small random perturbations of the equi- 
librium profile. The separation distance was increased 
in steps of s, taking the energy minimizing pattern of 
the previous step as the initial state in the CG scheme. 
To uncover the range of possible metastable pathways, 
we varied the step size s and, in addition, have con- 
sidered cases where the energy minimizing Fourier co- 
efficients (at do) arc perturbed randomly before being 
taken as initial choices for the next step. The per- 
turbations arc introduced by multiplying each Fourier 
coefficient with (l + r) where r is a random number be- 
tween — e and e, where e is called the noise amplitude. 

Fig. O'V depicts typical changes in the film morphol- 
ogy during the process of pull off starting from the 
critical distance dc where the instability originates. 
This fully nonlinear simulation (without any imposed 
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noise) shows the columns/cavities being laterally sep- 
arated by ~ 3/i at all separation distances until the 
maximum force is reached. Thus, the initial contact 
{d < dc) nanocavities predicted by the linear theory 
persist during puUofT. Figs.|2|3 and|2P show, respec- 
tively, the variation of contact area [a = area of con- 
tacts/total film surface area) and the force F on the 
contactor, with gap thickness for different step sizes, 
O.lnm, l.Onm and 2.nm. For small step sizes (< O.lnm) 
the debonding pathway is such that the configuration 
is trapped in the energy minimum corresponding to the 
initial instability. The contact area remains constant 
and the force increases almost linearly until a catas- 
trophic snap-off of the bridging columns. Remarkably, 
the maximum force Fmax that can be sustained before 
debonding is about an order of magnitude smaller than 
the maximum adhesive force Fj^^*. Clearly, debond- 
ing does not occur by a uniform detachment of the 
contacts, but rather by a different pathway, requiring 
much smaller pull off force, made possible by the initial 
pattern formation. The formation of bridges and cav- 
ities allows very high concentration of elastic stresses 
near the edges of the columns. For small step sizes, 
the elastic stresses build up to very high levels com- 
parable to the maximum adhesive force, without any 
intermediate small relaxations, since the structure is 
trapped in the original deep energy minimum. This en- 
genders a catastrophic adhesive failure for small step 
sizes. In contrast, larger step sizes force the struc- 
ture to hop through a succession of metastable states 
with lower barrier heights releasing energy intermit- 
tently leading to a continuous decrease in the contact 
area. The stresses at the edges of the contacts are not 
large enough to cause catastrophic detachment, but are 
sufficient to sustain peeling. For intermediate step sizes 
(curve (2) of fig. 133), the initial phase of pulling repro- 
duces the features of small step size results, followed by 
the large step size behavior. The escape from the initial 
high barrier state occurs only after after some stretch- 
ing of columns leading the release of pent-up elastic 
energy. The ascending branch ("elastic branch") of 
the force curve Fig. 3(C) reflects the linear increase 
of elastic stresses in the columns without any change 
in the contact area. The initiation of peeling limits 
the maximum force, after which it declines ("release 
branch") with further increase in the separation dis- 
tance and a concurrent reduction in the contact area. 
The release branch of the force curve is realizable only 
in displacement controlled experiments. 

A particularly simple model, which shows the essen- 
tial physics of the linear decrease of the area, approxi- 
mates the total energy as ^^iu'^{a/ {l—a)h)+aU (dQ—u) 
where a is the fractional contact area. For a given 
do I the minimum of energy occurs when the fractional 
contact area is a{do) « 1 - {3^J./2h\U{de)\)^/'^ do- This 
linear decrease in the contact area with increased sep- 
aration shows that debonding, even by the application 




FIG. 4: (A) Noise-induced, column-splitting pathway of 
adhesive failure for a more strongly adherent film (same 
as fig. H except, A = 10"^^ J and e = 0.01) (B) Variation 
of force F/fi with separation distances for various levels of 
noise (e = 0.0,0.001,0.01). 
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FIG. 5: Variation of snap-off distance with non-dimensional 
stiffness ratio. The slope of each curve is represented by p. 

of a purely normal force, actually proceeds by peeling 
of the contacts. Peeling from the contact edges requires 
much smaller energy penalty (and force) as compared 
to homogeneous debonding of flat contact areas. 

Although the peeling mode remains the dominant 
mode of debonding for small levels of noise, another 
pathway of debonding in the form of cavitation within 
the contact area leading to column splitting also ap- 
pears for high noise amphtudes (fig.^^), i. e., starting 
from initial conditions that are far from the solution 
branch being followed in the absence of noise. The 
column splitting mode is favored for higher adhesive 
strength and for more compliant films (higher value 
of AGh/fj,) where even smaller amounts of noise can 
induce this transition. Column splitting, when it oc- 
curs, results in precipitous decrease in the force, usu- 
ally followed by a regime of more nearly constant force 
(fig-Ej3)- Continuous peeling from the sides of the split 
columns prevents the build-up of elastic force in the 
constant force regime. This helps understand experi- 
mental observations of constant force regime which be- 
comes more prominent on rough surfaces that allow 
cavity initiation within the contact zones 0, 0| . 

The maximum pull-off distance was obtained for a 



wide range of parameters ft, ~ 0.1 — 50/im, ^ ~ 10 — 
10-21 J and A< ~ 0.1 - lOMPa, step sizes (s ~ 0.1 - 
2.0)nm, noise amplitude (e ~ 0.001 - 0.01) and AG - 
1 — lOOmJ/m^. Interestingly, the dependence of the 
puUofF distance dp on AG, ji and h is represented by a 
master curve of the form (fig. [S]) 
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where the nondimensional parameter {AG/d1)/{ii/h) 
is the ratio of the stiffness of the interaction potential 
and the elastic stiffness of the film. The exponent p is 
close to 1 for noiseless cases with small step size and 
decreases with increasing step size as well as increasing 
level of noise (the minimum exponent is about 0.6). 
As argued earlier, increased step size and noise levels 
can induce debonding at smaller distance by cascading 
through higher energy metastable states. 

Further, the maximum force from simulations was 
found to scale as 
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where exponents 7, S are —0.63 and 0.24, respectively 
for noiseless case with small step size (s = O.lnm), 
and the prefactor C = 38.5. In this case, F^ax ^ 
(AG)°-25/i-o.4^o.75^ g^lY other cases of higher step 
size and noise considered, 7 is found to be remark- 
ably constant at —1 ± 0.01 and the exponent S is 
0.8 ± 0.18. The prefactor for these exponents is in 
the range of 0.01 to 0.1. The force in this case scales 
as Fjnax ^ (AG)'''^/i-°'2/i'''2 showing a stronger de- 
pendence on adhesive energy but a weaker dependence 
on film thickness and shear modulus compared to the 
noiseless and small step size case. The above consider- 
ations (fig. ISJ also explain the long debated contention 
that the surface energy of soft solids as measured from 
debonding experiments is a non-equilibrium and non- 
unique property. 

The values of the maximum force per unit area re- 
quired for debonding are much smaller than predicted 
for debonding for flat surfaces {F/J^* ~ AG/de). For 
example with A = 10"^" J, h = 10.0/im, /i = O.bMPa), 
pf^f^ = 80MPa, Fynax = O.lMPa for step size Inm 
without noise and Fmax = 0.04MPa with noise. The 
ratio, F^ax/Fil^l ~ {pidl/ AGhY{h/ d,)^ , [n > m), 
where n is close to 0.2 and m close to for cases with 
noise and relatively large steps. This shows the dis- 
crepancy in the forces between the fiat and instability 
controlled modes of failure increases with decreasing 
shear modulus, increasing adhesive strength and film 
thickness. 

Representative 2D simulations (such as shown in 
fig.EJ) also confirm the underlying physics and the other 
results reported here. A detailed account will be pub- 
lished elsewhere. 




FIG. 6: The first two frames represent the instability pat- 
tern during approach (do = 5.72nm and 1.43nm) and 
the last two represent pattern evolution during withdrawal 
(do = 35.43nm and 43.43nm) The darker regions are the 
contact zones and the lighter shades denote the cavities 
formed within the film. During the approach a labyrinth 
pattern is transformed into isolated nanocavities which 
upon pulloff grow and fragment the contact zones. The 
wavelength of the pattern (~ 3h) also agrees with linear 
theory and experimental observations 



This letter resolves important open questions regard- 
ing the mechanisms and pathways of debonding at soft 
interfaces. The main results include, (a) the physical 
origins of the adhesion-debonding hysteresis, (b) for- 
mation and persistence of regularly arranged cavities 
and bridges during debonding, (c) metastable path- 
ways of debonding such as column collapse, column 
peeling and column splitting that require much larger 
pulloff distances and much smaller debonding forces as 
compared to detachment of flat surfaces, (d) complete 
quantitative dependance of pull-off distance and force 
on adhesion energy, shear modulus and film thickness. 
Formation of cavities engenders extremely high stresses 
near the column edges leading to the peeling of contact 
zones at much smaller average stresses than the adhe- 
sive strength. This is analogous to defects (dislocations 
and cracks) in solids which give rise to observed yield 
stress and strength much smaller than ideal values. 
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